Lattice QCD at non-zero isospin chemical potential 
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Abstract. Systems of non-zero isospin chemical potential are studied from a canonical 
approach by computing correlation functions with the quantum numbers of N ir + 's (Cnw)- 
In order to reduce the number of contractions required in calculating Cntv for a large N in the 
Wick's theorem, we constructed a few new algorithms. With these new algorithms, systems with 
isospin charge up to 72 are investigated on three anisotropic gauge ensembles with a pion mass 
of 390 MeV, and with lattice spatial extents L ~ 2.0, 2.5, 3.0 fm. The largest isospin density 
of pi ~ 9 fm -3 is achieved in the smallest volume, and the QCD phase diagram is investigated 
at a fixed low temperature at varying isospin chemical potentials, < fii < 4.5 m n . By 
investigating the behaviour of the extracted energy density of the system at different isospin 
chemical potentials, we numerically identified the conjectured transition to a Bose-Einstein 
condensation state at > m^. 



1. Introduction 

An important goal of nuclear physics is to investigate interactions between hadrons; however 
studying many body systems are hampered by the factorially growing number of Wick 
contractions naively required in computing the corresponding correlation function. In order to 
study multi- meson systems, a few algorithms [TJ [2J have been constructed to alleviate the cost of 
computing larger number of independent contractions. With these new algorithms, many-pion 
systems have been studied in Ref . [21 H] , and many-kaon systems have been studied in Ref . [5] , 
furthermore systems of mixed species have also been studied in Ref. [6] . In order to study more 
complicated multi-baryon systems, various methods have also been introduced. With methods 
constructed in [7], correlation functions of 28 Si have been computed within a manageable amount 
of time. 

Knowledge of systems at high density, with non-zero chemical potential, and non-zero 
temperature are vital to explore the QCD phase diagram. Systems with zero chemical potential 
have been investigated for a wide range of temperatures, and the transition from confined 
state at low temperature to the de-confined state at high temperature has been investigated 
substantially [8j. However, direct simulations of non-zero chemical potential systems are much 
harder due to the famous sign problem resulting from the non positive fermion determinants. 
But the sign problem does not exist for isospin chemical potential systems. Although QCD 
phase diagrams of these two systems are not the same, they do share some common properties, 
and it is also theoretical interesting to study non-zero isospin chemical potential systems. 

In this proceeding, systems of non-zero isospin chemical potential are studied from a canonical 
approach by directly computing correlation functions of n-7r + systems. By appling algorithms 
constructed in [2J, we studied systems containing up to 72 7r + 's on three anisotropic 2 + 1 flavor 



dynamical gauge ensembles, 16 3 x 128, 20 3 x 256 and 24 3 x 128, with £ = a s /at = 3.5, a s = 0.125 
fm, and with a poin mass of m n = 390 MeV and a kaon mass of mx = 540 MeV. The QCD phase 
diagram at a fixed low temperature, T ~ 20 MeV, for a range of isospin chemical potentials, 
m-w < Ml < 4.5 77% , are investigated, and evidence for a phase transition from pion gasses to 
Bose-Einsetin Condensate (BEC) state around fij = 1.3m n is numerically identified. 

The layout of this proceeding is as follows. In Sec{2] we briefly discuss the method applied 
to study n-n systems. Main results from the simulation are presented in Secj3j and we conclude 
in Secgl 



2. Methodology 

Non-zero isospin chemical systems can be investigated from the study of n-ir systems by directly 
computing their correlation functions. Because of the Pauli principle, the largest number of 
pions can be studied from a single source location is N S N C = 12. In order to study systems of 
more than 12 it's, additional source locations are required. The correlation function for a system 
of total n = Yli=i n i 7 r + 's with rij 7r + 's from source locations (y,, 0) is defined as: 

C ni ,...,n N (t) = ( ( J2 7V+(x,t)) (n-(yi,0)) . . . (V(y N , 0) ) ), (1) 



where the interpolating operator 7r + (x, t) = d(x, i)75«(x, t) and 7r~(x, t) = u(x, i)7sd(x, t). 

Computing the above correlation function includes 0(nl) contractions. For larger n and N, 
number of ways to distribute pions over N sources grows exponentially, which further increase 
the computational cost of all C ni) ... ntf 'a. However, correlation functions of n pions, C— , without 
knowing details of the distributions of ni can be constructed, which significantly reduces the 
number of correlation functions to compute as energies extracted from Cn is same as those 
extracted from C ni) ... jnjv 's as long as n = J2iLi n i- The explicit methods constructed in Ref.[2] 
are used to calculate Cn in this proceeding. Detailed discussions of the method and evaluations 
of n-ir correlation functions in momentum space can be found in Ref . [5] • 



3. Lattice Results 

Computations are performed on ensembles of nj = 2 + 1 anisotropic gauge configurations with 
clover-improved fermions generated by the Hadron Spectrum Collaboration and the Nuclear 
physics with Lattice QCD collaboration. Pion correlations function, C n7T (t), for n = 1,2, ■ ■ ■ , 72, 
are computed from three lattice ensembles of L 3 x T = 16 3 x 128, 20 3 x 256, 24 3 x 128. More 
details of these gauge ensembles can be found in Ref. [9] . 

The expected form of C n7r (t) computed on a lattice with temporal extent T is [6] 

LfJ / \ 

Cn.it) = ^ i n )A^Z^ E "-- +E ^ T / 2 co S h((E n . m -E m )(t-T/2)) + ..., (2) 
^ W 

where = 1 when m = n/2, otherwise = 2. E m is the ground state energy of a m-n 
system, the are overlap factors of thermal states with m tt's propagating backward around 
the temporal boundary [[] and the ellipsis denotes contributions from excited states. For infinite 
T, all thermal states vanish, and only the ground state survives, however the ground state 
is always contaminated by thermal states for finite T. Even for Ci2n computed on T = 128 
ensembles, the ground state does not dominate in any region. Under the current precision, C nn 
computed from ensembles with a temporal extent T = 256 shows a clear plateau region even 
for n = 72. Thus, in order to reduce thermal contaminations, the anti-periodic plus periodic 

1 When m = 0, Zq denotes the ground state contribution. 




Figure 1. On the left plot, ground state energies of n-Tr(E nn ) systems extracted from ensembles 
16 3 x 128 (red), 24 3 x 128 (blue) and 20 3 x 256 ( green) are shown. The black line represents the 
free energy of n non-interacting pions. Energy densities, e(pi), computed from three ensemble 
are compared on the right plot. 



propagation (A ± P) method [2] has been applied to effectively double the temporal extent of 
T = 128 ensembles. 

3.1. Ground state energies 

Ground state energies of n-7r systems are extracted by fitting a single exponential to C n7r (t) 
within time slices where a clear plateau can be identified in the corresponding effective mass plot. 
Statistical uncertainties are evaluated from the bootstrap method, and systematic uncertainties 
are from averaging systematic uncertainties computed on each bootstrap sample by moving the 
fitting window forward and backward two time slices. Extracted ground state energies from three 
volumes are shown in the left plot of Fig. [T] and the corresponding energy densities, e = E/V, 
are compared in the right. Because of the repulsive interactions between pions, E nn is larger 
in smaller volumes as pions are closer to each other, however energy densities computed from 
three volumes are approximately the same. 

3.2. Isospin chemical potential 

From the ground state energies of n-n systems (E nn ), an effective isospin chemical potential 
fj>i(n) = ^ can be computed from a finite backward derivative, m(n) = E n7T — ■E( n _iw. In 
order to take into account the correlation between E nn and E^ n _^ w extracted from the same 
bootstrap sample, ni(n) is also evaluated on each bootstrap sample. The extracted isospin 
chemical potential as a function of isospin density is plotted in Fig. [2] 

At small isospin density, m grows in an accelerating rate, agreeing with the prediction from 
X-PT [10], however it starts to deviate from the x~PT and begins to grow in a decelerating rate 
at large isospin densities, and at even larger isospin density it begins to flatten out. 

3.3. QCD phase diagram 

The change of the behavior of the isospin chemical potential at different isospin density signals 
a possible change of the physical state of the n-ir system. In order to investigate this change, we 
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Figure 2. The isospin chemical potential, ///, is plotted as a function of the isospin density, 
pi, from three lattice ensembles, 16 3 x 128 (red, pi = [0,9]), 20 3 x 256 (green, pi = [0,4.7]) and 
24 3 x 128 (Blue, pi = [0, 2.8]). The solid black line is from x-PT [TO] 



have studied the ratio of the energy density, e, to its zero temperature Stefan-Boltzmann limit, 



NfN r 

1^ 



(3) 



where Nf = 3 and N c = 3 is used in this formula. The ratio is plotted in Fig. [3j where it grows 
at small pi, reaches a peak around 1.3 m^, starts to decrease after the peak, and begins to 
flatten out at large isospin chemical potential. With a linear extrapolation of the peak position 
for different volumes as a function of the inverse of volume, the peak position in the infinite 
volume limit is Pp eak = 1.30(7) m^. 

At low temperature, as soon as pi reaches m n , one pion is produced out of vacuum, and the 
system is an interacting pion gas. When pi reaches the peak at pi = 1.30(7) m n , pions start to 
condense and the n-ir system becomes Bose-Einsten Condenstate (BEC). Such a transition to 
BEC state at pi > m n are also observed in two flavor QCD [11], and similar behavior occurs in 
two-color QCD [12]. 



4. Conclusion 

In this proceeding, we studied n-n systems from a canonical approach by explicitly computing 
correlation function of the n-ir system, C n (t). Ground state energies of the n-ir systems were 
extracted, and subsequently, the isospin chemical potential was computed as a function of the 
isospin density. 
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Figure 3. The e/esB is plotted as a function of fj>i/m n . 



By studying the ratio of the isospin energy density to its zero temperature Stefan-Bolzeman 
limit, the QCD phase diagram is investigated at a fixed low temperature, T = 20 MeV for a 
range of isospin chemical potentials, from fij = m n to /// = 4.5 m n . An evidence of transition 
from a pion gas to BEC is identified at /ij = 1.30(7) m n . 
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